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INTRODUCTION 

The chemical and physical Interactions causing molecular orientation 
at liquid crystal-solid Interfaces are a subject of fundamental physical 
Interest. They are also Important for practical applications in liquid 
crystal display devices. The successful operation of most electrooptic 
liquid crystal displays requires a well defined alignment of the liquid 
crystal molecules at the substrate surface. Special chemical or mechan- 
ical treatments of the substrate surface are required to obtain the 
alignment (Ref. 1). 

The optical behavior of an electrooptic liquid crystal display Is 
controlled by an applied electric field. A typical display geometry con- 
sists of a thin nematic layer (10-20 pm thick) sandwiched between tvo 
parallel surfaces, which have been treated to produce alignment. In the 
absence of electric (or magnetic) fields, the equilibrium bulk orienta- 
tion is determined by the elastic response of the liquid crystal to the 
interaction forces at the surface boundaries. The orientation of a 
nematic liquid crystal is conventionally described by the director, n, a 
unit vector representing the average orientation of the long molecular 
axes in a local region of space whose dimensions are small compared to 
the long range of the ordering forces present in a liquid crystal. An 
ordered nematic liquid crystal behaves optically like a uniaxial single 
crystal with the optic axis parallel to the director tt. For liquid 
crystals with positive dielectric anisotropy, the application of an 
electric field tends to reorient the molecules along the field direction. 


The onset of reorientation occurs at a critical voltage, producing a 
sharp transition in the optical properties. 

Among the electrooptical qualities desired for good device perform- 
ance are sharpness of transition and transparency. These qualities are 
critically Influenced by the direction and uniformity of the initial 
surface alignment and by rHe strength of interaction (pinning strength) 
between the surface forces and the liquid crystal molecules (Ref. 2). A 
measure of the pinning strength is the anchoring energy. Quantitative 
evaluation of the anchoring energies for various surface treatments, 
types of alignment and different liquid crystal materials is needed for 
a better understanding of the diverse factors contributing to surface 
alignment and could be useful to the development of display technology. 

The most commonly used model for the surface energy (Ref. 2), re- 
flecting the anisotropic nature of the surface interactions, has the form 
0 ) ■ W sin^(A0) (1) 

where m (ergs/cra^) is the energy required to rotate the director away 
from its preferred surface orientation (called the "easy direction") 
through a small angle A0. The coefficient W in this expression is the 
surface anchoring energy. Equation (1) has been used in measurements of 
surface energies corresponding' to rotations of the director both in a 
vertical direction away from the plane of the substrate and in the plane 
of the substrate. Extremes of observed anchoring energies, corresponding 
to different surface treatments, range from about 10"^ergs/cm^ (weak 
anchoring) for rotations in a vertical direction up to about one erg/cm^ 
(strong anchoiring) for rotations in the plane of the substrate (Ref. 2). 


Another parameter sometimes used to characterize surface energies is 
"extrapolation length" (Ref. 3, p. 74), whose magnitude is on the order 
of K/W, where K is an elastic constant of the liquid crystal (gener- 
ally about 10“^ dynes). The extrapolation length defines an effective 

sample thickness when a distortion has been Imposed on the bulk orlenta- 

// 

tlon and can be regarded as a measure- of the continuation into the sub- 
strate of the angular variation of the director. For strong anchoring 
conditions, in which the surface energies are comparable to or larger 
than the interaction energies between molecules, the extrapolation length 
is on the order of or smaller than average molecular dimensions (.^100 A) 
and for all practical purposes Is essentially zero, For weak anchoring 
conditions, the surface energies are much smaller than the Intermolecular 
interaction energies. Then the extrapolation length is much larger than 
the molecular dimensions and can be as high as 100 pm. 

Only a few methods exist for the experimental determination of sur- 
face anchoring energies. In one of these methods (Ref. 4), quantitative 
estimates of anchoring energies have been obtained from a combination of 
optical analyses of the variation in director orientation a^iiross surface 
dlscllnatlon lines and measurements of the equilibrium line widths. Sur- 
face dlscllnation lines are a common defect observed in liquid crystals 
and are lines of discontinuity in the molecular orientation attached to 
the substrate surface. In a second method, analysis of the wall effects 
on the magnetic Preederickse transition in a homeotropically aligned 
nematic cell have been used to evaluate anchoring energies for substrates 
treated with various surfactants (Ref. 5). In yet another method, 
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anchoring energies were determined in twisted nematic cells from measure- 
ments of the director rotation caused by application of a raagnetlc field 
(Ref. 6). 


A new method for evaluating surface anchoring energies from measure- 
ments of the director orientation at the surface boundary of a nematic 
liquid crystal cell has been presented recently by Riviere, Levy and 
Guyon (Ref. 7). The theoretical and experimental results of Riviere, 

Levy and Guyon (RLG) are briefly reviewed below. 

The sample used in the RLG experiment consisted of a slightly wedge- 
shaped nematic liquid crystal cell formed by two nearly parallel glass 
plates. The sample geometry is shown in Figure 1. 



Figure 1. Liquid crystal sample geometry used for determination 
of surface anchoring energy from tilt angle measurement. 

The angles (|)j^ and <|>2 are the easy directions of alignment and repre- 
sent What the tilt angles would be for each surface alone. The angles 
and 02 actual tilt angles and are distortions induced by the 

opposite surface. The small wedge angle a in the RLG study was 18 
minutes of arc. The upper surface was treated with an organosilane to 
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produce a strong homeotroplc alignment (i >2 ~ 90°) . The lower surface was 
coated with an obliquely deposited SiO film which provided an easy dlrec'* 
tlon of ~ 30°. The angle 0 was the angle formed by the director n 
and the X-axis. Rotations of the director were assumed confined to the 
XOZ plane. 

To derive an expression relating surface anchoring energy to the 

<■/ 

tilt angle and sample thickness d, RLG started with the free energy 
density as given by the isotropic elastic continuum theory for liquid 
crystals, 

Fd * I K (76)^ (Z) 


To simplify calculations, the one-constant approximation was assumed, in 
which the liquid crystal is considered as an isotropic elastic medium 
whose elastic properties are represented by a single elastic constant, K. 
For this free energy density, the Euler-Lagrange equation for minimisa- 
tion of the total bulk free energy is simply Laplace's equation, V^0 = 0. 
Because of the small wedge angle a, it was also assumed that the varia- 
tion of the director angle 0 in the X-direction was much smaller than 


in the Z-direction, viz., ■!“ << 


For this assumption V '0 


3^0 

3z2 


Integration of the one dimensional Laplace equation gives 
6 = az + 0„ 


0 . 


(3) 


30 


" ®1 


where d is, the sample 


where a = For a given x, a - ^ 

thickness at that point, The surface energy at the lower surface 
(subscript 1) was assumed to be given by 


“SI = sin (01 - W 


(4) 
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SubsCititloa of Eq. (4) Into the equation representing the balance 
between the surface and elastic torques, 


‘"“si 

de 



(5) 


resulted in the following expression for the surface anchoring energy 

K (82 - 8l) 

“ d ain[2(6i - ti)] 


To enable evaluation of from experimental measurements of the tilt 

angle 0j^ and sample thickness d only, some additional approximations 

were made. One approximation was that ^ = M and that 

d d d 

&(p was constant over the region of measurement. Also for small angles 
(®1 ■ 4’l)» it was assumed that sin[2(0j^ ^ ^ 2(0x - ^i)« Using tilt 

angles, 0| and 6*[, and the corresponding sample thicknesses, d’ and 
d*V, measured at two Xoeatlons along the lower surface, Eq. (6) can be 
rewritten as 


”i “ 2(0*^ - ep \d' ■ 

The tilt angles were deduced to within ±0.3° from reflectivity measure- 
ments at the critical angle of parallel polarized light from a laser 

» t 
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light source (Refs. 7»8). Using an approximate value of K'^lO erg-cm, 
RLG reported a calculated value for of (2.1 + 0.8)xl0 ergs/cm , 
which generally agreed with previously reported values of the surface 

I 

anchoring energy for the same surface treatment. 

The purpose of this thesis is Lo extend and generalize the theory 
presented by RT/! to Ittclude non-eciiial elastic constants and a two 



dimensional variation of the director orientation. The sample geometry 
and parameters are the same as defined by Figure 1. The sequence of 
presentation is as follows. First, there will be a very brief discussion 
of relevant aspects of the elastic continuum theory for nematic liquid 
crystals. Then, an analytical solution to the one-dimensional problem 
allowing for elastic anisotropy, i.e., assuming non-equal elastic con- 
stants, will be presented and discussed. This will be followed by a 
special solution for the two-dimensional problem, appropriate for the 
slightly perturbed one-dlmensional geometry representing the wedge- 
shaped liquid crystal sample. Next, a comparison will be made between 
the surface anchoring energies obtained using the one-dimensional, 
elastically isotropic approximation and the two solutions for elastic 
anisotropy presented here. Finally, an experiment for measuring surface 
anchoring energies using the analytic solution for the one-dimensional 
case will be proposed. 

FUSTIC continuum THEORY 

A number of large scale phenomena Involving the response of the bulk 
liquid crystal to external disturbances can be successfully described by 
the elastic continuum theory (Ref. 3, Chap. 3; Ref. 9, Chap. 8), in which 
the liquid crystal is treated as a continuous elastic medium. In the 
absence of applied electric or magnetic fields, the thermodynamic equilib- 
rium state is determined solely by the elastic response of the bulk liquid 
crystal to the surface interaction forces at the walls of its container. 
The Surface forces induce a static distortion of the director orientation. 
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The starting point for calculations is the equation for the elsstle 
free energy density, expressed in terms of the director field, n(r). The 
elastic free energy density is given by 


k 


V'n(r) 




n(r)’V X 


n(r) 


+ K. 


A A . 

n(r) X (V X n(r)) 


1 


( 8 ) 


Here Kj^, K 2 and K 3 are the Pranck elastic constants and correspond to 
the three basic types of distortion for a nematic liquid crystal, namely, 
splay, twist and bend, respectively. The basic problem in applications 
of the elastic continuum theory is to determine that director configura- 
tion which minimizes the total free energy, The total free enargy is 
given by 

F-Fg + Fg 

“ ^ (i =^1,2,3) (9) 

Volume 
of sample 

where Fg is the total bulk elastic free energy and Fg is the total 
surface energy. The equilibrium conditions are determined by minimizing 
the total free energy with respect to all variations of the director n(r) 
subject to the constraint n*n - 1. In actual calculations it is more 

A 

convenient to write n in terms of appropriate polar angles and minimize 
F with respect to all variations of these angles# 

For the sample geometry shown in Figure 1, all rotations of the 
director are assumed to take place only in the XOZ plane. There is then 
only one polar angle, 9, as defined in Figure 1. In terms of 0, the 

A 

components of n are 




« n’Cjj • COP 0 

Ity “ ft’Gy « 0 

nj 5 '« n'% « sin O 


where e^, Oy» vectors for the Cartesliin coordinate 

systemi By direct calculation 
V'n " - sin 00JJ + cos 00 ^ 
and 

Vxn “ - (cos 00JJ + sin 005 ,)e 
where 


£. 30 o 30 

«x - 9-x> H 




Then 


(V*n)^ * (sin 00j^ - cos 00^)^ 
n*7 X n ® 0 


Cn 


n X (V X tt) 


% 2 
“ (cos 00„ + sin 00 „) 

-X. « 


The free energy density Is 

'd ' \ {h^I + »2Vx + %®x) 

where 

gj^ = K2^sin20 + K3COS20 
g2 == 2(K]^ - K3)sin 0 cos 0 
S3 = K^cos^e + K3sin26 


(10) 


( 11 ) 
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In the one-constant approximation whora K Kj^ *• %, Bq. (10) roducea Co 
the F^j given by Eq. (2)^ 

Using Eq* (10) for che elastle free energy density, Eq* (9) for che 
total free energy eon be written 


F w 


/// 


,0x>0a)dxdydz + / / 

i»l , 2 



(12) 


where 


“l W^sln^(0^ - «|>jl) |: 

and doj^ Is an elemenf‘\of area on the ith surface. 
The minltnization of F requires that 


6F 


/// 


2)dxdydz + 


D// 


6wida^ “ 0 


1 * 1,2 


(13) 


(14) 


for arbitrary variations <50 in the director angle 0. According to the 
calculus of variations, the variation (SP^j can be written as 

0Pd 0Fd 3Fd 
«Fj(e,9,,8^) «e + 30®- Mx + 8?: 

X z 


0F, 


A ae + i£d JL (50) + Jl, (66) 

ft Art A^_ ' Art A_ ' ^ 


3F. 


00 00JJ . 0x ' 302 


Ha.± 

30 3x 




+ X 60 1 

3z U02 


( 15 ) 


where (^0) and 60^ » ~ (50). 

The variation in the total Tree energy is then 



DeGennea has shown (Ref. 3, pp. 73-76) that when the surface forces 
are strong enough to orient the director n in a well-defined direction 
("easy direction") at the surface, the extrapolation length is on the 
order of molecular dimensions. It then suffices to neglect the surface 
terms and to minimize the bulk free energy term only, with fixed boundary 
conditions for n. This is a common practice in most calculations using 
the elastic continuum theory. In that case one would consider only the 
first integral in Eq. (16). As a condition for local bulk equilibrium, 
one has the familiar Euler-Lagrange equation, 

la n7i 

36 ■ 3x ^38x7 ” 3“ ° 

In the case considered here, however, surface energy plays an impor- 

tant role. We are ;lot5king( at the influence of one surface, manifesting 

'W f . . 

Itself in torques tiansmitted through the liquid crystal, in modifying 
the orientation of the director at the opi>osite surface. The surface 
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terms cannot be Ignored, In addition to Bq. (17), which expreaseg the 
condition for local equilibrium In the bulk> we hove the following con- 
dition for mechanical equilibrium at the substrute surfocoB, 


X^X y 


dfe.| ** 0 


(18) 


As will be seen, in subsequent applications to the one- and two- 
dimensional cases, Bq, (18) introduces in a natural way the boundary 
conditions expressing the balance between surface energy torques 

36 0«0 


exerted on the liquid crystal and the elastic torques exerted 


by the liquid crystal. 


VARIATION IN ONE DIMENSION 

We consider first the case where it is assumed that variations of 


the director angle 0 in the X-direction can be neglected in the calcu- 
lation of the equilibrium conditions, i.e, , « Q^. This assumption is 
equivalent to taking the two surfaces to be parallel. Equation (10) for 
the free energy density is then 


F,i{3.0,) - A 836^ 


(19) 


2 2 

where gg ^ *lino® 6 K 3 sin 0. The Euler-Lagrange equation for bulk 
equilibrium is 



( 20 ) 
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Substitution of Eq. (19) Into Eq. (20) gives for the Euler-tagrange 
equation » 


2 836^^ - 0 


( 21 ) 


where and 9 “ With z being the only Independent 

“39 90 zz 

variable, 9^ and are total derivatives. Eq. (21) is an ordinary 

second order differential equation which can be integrated directly as 
follows. First multiply by Sg. Then 


2 ®39®z ®3®z® 


zz dz V 2 


^ (I ^A) 


( 22 ) 


dz 


(Fd) 


* 0 

The free energy density is thus a constant which can be evaluated in 
terms of the boundary conditions as follows. Let Fjj - C. Then 

1/2 


d9 

dz 


■(«) 


Integration of Eq. (23) gives 



flic (32 - z^) 


Let d = Z£ "• z^ be the sample thickness and 


(23) 


(24) 
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El(8l. 


gi«de 


/®2 

<h 


cos^e -f K2Sin^6)^^^d0 


Ei(6j^>62) Is an elliptic type Integral and Is readily evaluated, given 
0 ]^, 02» Kj^ and K3. The subscript l*in‘ '£3(02,02) refers to the one- 
dlmenslonal variation of 0 . The free energy density Is then 




For two parallel surfaces, the elastic free energy density given by \ 

" • . . t 

Eq. ( 26 ) is constant throughout the sample. For the slightly wedge*- 
shaped sample, calculation of the equlllblrum conditions based on the 
assumption that 0^ can be neglected indicates that F^ remains 
constant (or approximately so) . This Implies that at a given location 
X, where the sample thickness is d, 0^ and 02 should assume values 
such that Eq. ( 26 ) Is satisfied. 

The surface terms In the variational calculation, Eq. ( 18 ), give the 
equilibrium conditions for the balance of torques at the surface. From 
the balance of torque equation, one then obtains an expression relating 
surface anchoring energy to measurable parameters such as sample thick- 
ness d and tilt angles 63 and 62. 

Consider the lower surface, z * Z3. The variation in surface energy, 
3 uJl j 

«“i • 1;- L ^ - 1= 
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■ 2Wj^sln(e - (Jij)cos(e - <|)j)6' 




Wj^sin 2(6 - (j)j^)6e 




(27) 


For the one-dlmenaional case, where — = « 0, the first integral in 

30X 


Eq» (18) is 
,z 
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(28) 


Combining the lower limit of Eq. (28) with the variation given by 

Eq. (27), one obtains the following equilibrium condition expressing the 
balance of torques at the lower surface, 


30, 


Z“Zi 


+ Wj^sin[2(e - (j)^^)] = 0 

I Z=ZJ^ 


(29) 


From Eqs. (19) and (26), 


Hi 

30, 




(30) 


The anchoring energy Wj^ at the ’lower surface is then 


W, 


1/2 

d sin[2(0j^ - 0j^)} 


1/2 


(Kj^cos^Qj^ + K3Sln20j^) (03^,02) 

d sinj2(0j^ - 


( 31 ) 


1/2 1/2 

In the one-constant approximation where 

E (01,02) “ K^'"^(02 - 0x>- is then identical to Eq. (6). An 

expression for the surface anchoring energy W 2 at the upper surface 
can be obtained In a similar manner. 

VARIATION IN TWO DIMENSIONS 


For variations of the director orientation angle 0 in both the x 
and z directions, the free energy density is given by Eq. (10) as 


Fd(0.0x»6z) = ^ (Sl^x + 820 x 02 + 830 z) 


( 10 ) 


where g^, g 2 and g^ are the functions of 0 only defined by Eq. (11)^ 
The total free energy, elastic and surface, is given by Eq. (12) as 


F - 


/// 


Fd(0,0x.0z>dxdydz + 




( 12 ) 


1 * 1,2 


The Euler-Lagrange equation representing the bulk equilibrium condition 
is given by Eq. (17) as 



Upon substitution of Eq. (10) into Eq. (17), the Euler-Lagrange equation 
becomes J' 


f (0,0x>0z*0xx’®xz»®zz^ - 2 (8i00x + 8200 xOz 

+ (810XX + 820X0Z + 83822) = 0 (32) 
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where 


^ (i » 1,2,3), 0XX * 


ae 


2’ xz 


2 

3x3 z ^ ®zz 


2 

3_ e 

. 2 * 


3x^ "" ““ 3z' 

Eq. (32) Is a "quasi-linear" second order partial differential equation, 
so called because of being linear In the highest order derivatives. A 
general solution to this equation was not found In the literature. 
Attempts to obtain such a solution by systematic or ad hoc methods were 
unsuccessful (Ref. 10). 

A special solution to Eq. (32), similar to the solution to the one- 
dimensional Euler-Lagrange equation, Eq. (20), can be derived. This 
"linearized" solution is appropriate for the slightly wedge-shaped sample 
geometry, which is only a very small departure from the one-dimensional 


parallel surface geometry. 

An insight into how to approach a solution can be gotten by looking 
back to the one-dimensional case. For the two surfaces parallel, the 
calculated free energy density, Fj, was found to be constant throughout 
the sample. This was an exact result. The same calculation was assumed 
approximately valid for the slightly wedge-shaped geometry with the 
implied result that the free energy density remained approximately con- 
stant over the dimensions of the sample. In the one-elastic constant 
approximation, which gives qualitatively the same results as Ki Ko, 

r •< X ^ 

l2 0o*“0i 

Fj “ K02 with 0z = — i. The difference in thickness between two 

points along the sample is Ad = tan a Ax = ^ Ax. For tan a « 1, 

Ad << Ax. Thus d varies very slowly with X and consequently So do 
02 and Fj. Therefore, for >> 1, the assumption that the free energy 
density is constant over the dimensions of the sample is a reasonably good 


approximation. 
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Although variations of 0 with both x and z will be allowed, 
let us assume that Fj remains eonstant over the volume of the sample. 
This assumption will enable a linear solution to Eq. (32). Let 

Fd(®*®x>®z) “ *2 (8l®x 82 ®x®z 83®z) “ ^ (33) 


where C is a constant to be evaluated in terms of the boundary 
conditions. 

For first order partial differential equations in which the inde- 
pendent variables do not appear explicitly, of which Eq. (33) is an 
example, Charpit’s method (Ref. 11, pp. 69-73) can often be effectively 
applied to obtain a solution. In this case, two of Charpit's equations 
assume the especially simple form 


d0. 


de, 


(34) 


where d0„ and d0„ are differentials. It follows immediately that 


®X 302 


(35) 


where a is a constant. Substituting Eq. (35) into Eq» (33), we find 
Fd(e, 0 x.e 2 ) ^ Fd(e, 02 ) and 

Fd(e.e,) - i (a^gj + aS2 + g3)e2 

= I GSj (36) 


= c 
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where 

G - + ag2 + (37) 

Then 



(38) 


The expression for in Eq. (38) is similar in form to the result for 

one-dimensional case given hy Eq. (23)- For a ■ 0, whence 0 * 0, the 
results are identical. It is easily verified that the solution to E(j. 
(36) also satisfies the Euler-Lagrange equation, Eq. (32). For example,! 
substitution of Eq. (38) into Eq. (32) shows for constant C that 


1 ^8l0®x ?20^x«z 830655) + (8165,^ + 820 x 5 . + 

“ f <^^«ie «820 + 83 o)(f) - (^hi + ag2 + 83)(^)(|f) 

« C 3G C 9G 
G 30 ~ G 30 

- 0 


A complete integral for 0 can be obtained by Integrating the total 
differential 
d0 * 0^dx + Qgdz 
* 0z(a dx + dz) 


/2c\1/2 

■ dx + dz) 

It is convenient at this point to define a new function 
s * a x +■ z 

Then ds = a dx + dz and 


(39) 

s(x,z) where 

(40) 
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The Introduction of the function s in Eq. (39) has in effect reduced 
the number of independent variables from two to one. To facilitate the 
integration of Eq. (42), let us look first at the geometrical signifi- 
cance of 8. For s ” constant, 

ds " a dx *f dz ^ 0 =$» fS. a -a 

dx 

Similarly for 6 constant, 
d6 = 6z(a dx + dz) = 0 gx 

In the XOZ plane, lines of constant s and 6 are parallel straight 
lines with slope -a. If one assumes rigid anchoring at the upper sur- 
face, i.e., 02 = <j >2 ® constant, the constant a is related directly to 
tan o, in fact, a - -tan a. Tljis can be directly verified in the follow- 
ing simple way. At the upper surface, 0 =02 constant implies that 
everywhere on the surface that V0 is normal to the surface, where 

A A 

V0 ■ O^e^ + 0z®z’ normal to the upper surface is given by 

A A 

N = sin a - cos a e^ 


Then 
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70 X N “ ^ ® 

"* ©^(a cos o + sin a)Oy 

- 0 

from which a •“ - tan a. The surfaces of constant 0 and s are thus 
planes parallel to the upper inclined surface. The gradients, which are 
normal to these surfaces, are related by 

A ft 90 

70 - 02 (a ex + e{.) « -g~ 7s (43) 

Generally, ©2 Is not expected to be constant over the dimensions of 
the sample, even for strong honieotroplc alignment at the upper surface, 
and most likely would vary slowly with x. For large differences between 
0j^ and © 2 * the parallel planes representing the surfaces of constant 0 
cannot intersect both the upper and lower surfaces. Consequently, 
ja( < tan o. In the special case of rigid anchoring at the upper surface 
just considered, ja| has a maximum value of tan o (5.24x10”^ for a “ 18 
minutes of arc). Generally, 0 < ja| < tan a. 

As in the one-dimensional case, the constant G denoting Che free 
energy density can be evaluated in terms of the boundary conditions. Let 
(xa*zia) he a point on the lower surface with tilt angle 0^^^ and 
be the opposite point on the upper surface with tilt angle 02g. Inte- 


grating Eq. (42) along ds = dz between these points, we obtain 





where dg Is t;he sample thickness at x " Xg. let 
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E2(0x.e2) 


/ 




1/2 

G ' do 


(45) 


Like Ei(Ox* 02) (25), E2(Oi,02) Is an elliptic integral and can be 

evaluated numerically, given a, Kj, K 2 , 0^ and subscript 2 

refers to the integral being a result of a two-dimensional analysis. 

With Che understanding that Oj^, O 2 and d all correspond to the same 
X, we drop the subscript a in Eq. (44). Combining Eq. (44) with 


Eq. (45), we obtain for the free energy density 
|2 


I|2(®1»®2)I1‘ 

2d2 


(46) 


For a ■ 0, Ej “ E 2 and Eq. (46) becomes identical to Eq. (26) for the 
one-dimensional case. 

An alternate approach to the solution of the Euler-Lagrange equa- 
tion, Eq. (32), which yields the same expression for the free energy 
density, Eq. (46), is the following. Because the slightly wedge-shaped 
sample geometry is such a small departure from the parallel surface 
geometry (tan a << 1), assume that 0 can be represented by a function 
of a linear combination of x ‘and z, ax + z, where a is on the order 
of tan a. Let 0 0(s) where 

s « ax + z (47) 

Then 


„ 80 9s 80 ^ 

* 8s 8z ’ "Ss’ ^ 


8J_ 8 s 
8s 8x 


80. 

8S 




( 48 ) 
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Substitution of Eq, (48) into Eq. (10) for the free energy density gives 


(». i) ■ i « 


where C is the function of 0 only defined by Eq. (37). Substitution 
of Eq. (48) into the Euler-bagrange equation, Eq. (32), gives 


1 3 /30\ 

prm y HMw l I 

2 3s\3s/ 


4- G - 0 

3s2 


which upon multiplication by ~ can be rewritten as 

3s 


2 3F. 


1 p /36V 

2 ^ W 


Let 0^ = 1 ^. Then 


dFd<0*6s) ' d0 + 


\30 3s 30g 3s/ 


Thus - G is constant and C 


r /iiV 


can be written in the form 


of Eq. (42) and integrated to obtain Eq. (46). 

As in the one-dlmenslonal case, the surface energy, can now be 
written in terms of the sample thickness, d, and the tilt angles 0^^ and 
02* The first step is to derive the equilibrium conditions for the 
balance of torques at the surface from Eq. (18). Consider the first 


Integral in Eq. (18), viz.. 
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( 7";^ Idxdyd! 

30x ) 92 \9055 /J 


Using 0x ** aOg and taking thn Indicated dorlvativmsi the Integrand 
becomes 


9 

fx 


(i^ i (IS' 


“ fz h CdO) - CC.0J,50) 


'2-' 32 
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The first integral in Eq. (18) is then 

/;[/: 


1^ I 


~ (G0j5fi0)d2 dxdy » 

'zi J 



[G02691 >dy 
21 


(51) 


(52) 


Using the lower limit in Eq. (52), the total variation of the surface 
terms at the lower surface is given by 

" [G0z6C|^^«xx‘^xdy 4* JJ ^ tDj^dxdy 

- fj (j}<50z + W] sin (2(0 - ♦J 00^ dxdy 

W W »- 2 “ 2 x 

- 0 (53) 

Eq* (53) , yields the equilibrium condition for the balance of torques at 
the lower surface, vis.. 


G0 


2 17, 


'21 


Wi8in[2(0i - 


, 1/2 


(54) 


Substituting G0z » (2CG)"'‘- and 2C 
surface energy at the lower surface 

" d sln[2(9x ■' ijix)] 


52 ( 01 > 02 ): I ' 


we obtain for the 


( 55 ) 
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which for a « 0 becomes Identical to the one“dimenslonol result given 

by Eq. (31). 

HEBULTS AND DISCUSSION 

Equations for the surface ancliorlng energy have been derived for two 
cases which allow for elastic anisotropy, l,e., i K 3 . These equations, 
Including the one-dlmenalonal Isotropic result of RLG, are summarized In 
Table 1 below. 

TABLE I. Surface anchoring energies, Wj^, for Isotropic and anisotropic 
elasticity. 




One-dlmenslon isotropic 
K, - Kj. “ K3 

K( 0 « - 0 ,) 

Eq. (6) 

d sin 1 2(01 - ©1^ 

One-dlmenslon anisotropic 
% i K3 

1/2 

(S3 )es0,El(0i,02) 

"■"■‘■"■■"""'P'"' — Eq. ( 31 ) 

d sin 12(03^ - (J>i)j 

Two-dlmenslon anisotropic 

Ki#K 3 


P Eq. (00) 

d Sin [ 2(01 - (^i)J 


The quantities $ 3 , E]_( 0 i, 02 )» G and E2(®1»®2) Table 1 are defined 
by Eqs. (11), (25), (37) and (45), respectively. 

For numerical comparison of elastic anisotropy effects on the calcu- 
lated values of Wj^, It suffices to use only the numerators of the 
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expressions given In Tnble I. Of the elostle constants, the bend con- 
stant Is the more sensitive Indicator of elastic anisotropy, particularly 
near the neraatlc-smectlc A transition. To empliaslse the dependence on 
Kj, the numerators ate rewritten as follows. 


1/2 1/2 

(83 )g-0j^%(®l*®2) * *^3% 

where £3 ■ 83 *• (1 + k cos^G) 


Ex<^1»02) “ E1(0x,03) 



and 



- 1 . 


Also 


(£^^^)e.9/2(01^92) - K3(Gl/2)g_Q^E2<0i.02) 


(56) 


(57) 


where G » G » 


(1 + a^) + k(cos 0 + a sin 0 )^ 


12(01,02) “ -jr ^2(01*02) “ J #^^d0 

^3 n.. 


/' 

•^01 


Equations (56) and (57), representing the one- and two-dimensional aniso- 
tropic cases, are to be compared with K(02 “ 0i) • Evaluation of as 
shovm requires knowledge of 0j^, 02 and d. The slightly wedge-shaped 
sample configuration used by RLG (Kef. 7) did not allow simultaneous 
accurate measurement of both 0i and 02- To enable evaluation of Wj|^ 
from measured values of 0i and d only, RLG made the assumption men- 
tioned in the Introduction, namely that (©2 ~ 0i) * (<i>2 ” "J*!) 6<p 
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remained approxlmaCely constant over the region measurement. Eq. ( 6 ) 
for then becomes 


KA(t> 

d sinfi( 0 i - 


(58) 


The above Assumption also predicted the experimentally observed approxi- 
mately linear variation of 9]^ with 1/d. This linear dependence is 
readily seen by rewriting Eq. (58) as follows, using 
sin[ 2 ( 6 x - <j>l)] = 2(01 - <{ji) for small angles ( 0 i - ^i) . 



For numerical evaluation of Eqs. (56) and (57), the assumption 


guiding the choice of 0 x and Og is that ©x vary linearly with 1 /d 
for small angles (0j^ - (jij^) . Equations (56) and (57) should then remain 
constant. The surface anchoring energy in the one-dimensional isotropic 
case (Ref. 7) was evaluated using ~ 30° and <fi 2 ~ 90° for 0^ and 

02 , respectively. The same procedure is followed to enable a numerical 
evaluation for the anisotropic cases presented here. The calculated 

i .] 

values of Eqs. (56) and (57) are presented in Table II for a range of 
anisotropic conditions, from the isotropic elastic case where K»K 2 ^ = K 3 
to large anisotropy where K 3 = 10 For comparison purposes, it was 

assumed K = Kj - 1, with K 3 allowed to vary. 

It was shown previously that |a| tan a. For the purposes of num- 

erical evaluation, it was assumed that a = -tan a. For ot equal to 18 
minutes of arc, tan a = 5.24xlO”3. Because of the smallness of this 
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factor, the two-dimensional results shown In Table 11 do not differ sig- 
nificantly from the one-dlmenslonal anisotropic results and hence will 
not be discussed further. 

As Table 11 shows, anisotropy effects are significant. Whereas, 
both £2 l.l(®l»®2) decrease with increasing K3, the overall effect 

resulting from the factor K3 Is an Increase in the calculated values of 
W3. The ratios of for the anisotropic case to W3 for the iso- 

tropic case range from 1.23 for K3 = 1.5 up to 4.83 for K3 = 10 Ki* 

, '} 

The above calculations show that inclusion of elastic anisotropy 
effects can contribute significantly to values of the surface anchoring 
energy calculated from measurements of surface tilt angles. In the ex- 
periment with the wedge-shaped sample (Ref. 7), lack of knowledge of 62 
hindered including elastic anisotropy effects directly in the calculation 
of surface energies, although one can assume an average of K3 and X\3 
for the value of K. The liquid crystal used in this experiment (Ref . J?) 
was 4-cyano-4 ' n-hexylbiphenyl (6CB) whose elastic constants have a fa tin 
of about K3/K1 = 2 at room temperature. Table II shows that at 
K.3 = 2Kj^, the anisotropy effect could cause almost a 50% difference in 
the calculated value of the anchoring energy. 

It would be desirable to rigorously include anisotropy effects in 
the calculation of Wj^. An experiment is now described which permits 
this and uses the exact one-dimensional equation for W^ (Eq. (31)). The 
sample configuration for this experiment consists of two parallel sur- 
faces where the nematic layer is of constant thickness, d, over the 
sample. The experimental technique of RLG (Refs. 7 and 8) is used to 
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tijeasure the director tilt angles (6j^ and ©2) et both surfaces. As 
djiscrlbed In the Introduction, the RLG technique consisted of deducing Oj^ 

from reflectivity measurements at the critical angle of parallel polar- 

% 

j! ized light from a laser light source. This teclinique is capable of de- 
i tectlng very small variations in liquid crystal orientation at the 
substrate surface. Tilt angles were measured to within ± 0.3 degrees, 

A glass prism of high index of refraction (about 1 . 9 ) was used for the 
lower substrate. In the experiment being proposed, use of a glass prism 
with appropriate index of refraction would enable measurement of 82 at 
the upper surface as well. Both 81 and 82 could then be measured at 
a number of locations over the sample and statistically averaged to 
obtain representative values for the surfaces. An advantage of the pro- 
posed experimental approach is that, by measuring both and 62, no 

approximations need be made in the calculation of the surface anchoring 
energy. Elastic anisotropy would be directly accounted for. The only 
assumption, which is necessary in any case, is the surface energy model 
in Eq. ( 1 ). 

The proposed experimental technique is potentially useful for the 
systematic determination of surface anchoring energies, strong and weak, 
for a variety of surface treatments and nematic materials. Such knowl- 
edge could contribute to both a better understanding of liquid crystal- 
substrate interactions and the development of liquid crystal display 


tecluiology. 


— !r ? 

i 

} 
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